Detection of gravitational wave bursts by interferometric detectors 
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We study in this paper some filters for the detection of burst-like signals in the data of interfero- 
metric gravitational-wave detectors. We present first two general (non-linear) filters with no a priori 
assumption on the waveforms to detect. A third filter, a peak correlator, is also introduced and 
permits to estimate the gain, when some prior information is known about the waveforms. We use 
the catalogue of supernova gravitational-wave signals built by Zwerger and Muller in order to have 
a benchmark of the performance of each filter and to compare to the performance of the optimal 
filter. The three filters could be a part of an on-line triggering in interferometric gravitational- wave 
detectors, specialised in the selection of burst events. 
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I. INTRODUCTION 



O 

Long baseline interferometric gravitational- wave detectors such as LIGO 0, VIRGO f|, GEO600 f| or TAMA300 
pj are now in their phase of construction, and should be fully operational in the first years of the next millenium. 
Sources of gravitational waves that are expected in the bandwidth of these detectors, all involve compact objects such 

in : 



as black holes (BH) or neutron stars (NS); see ^ or || for a review. Among them, inspiraling binaries seem to be the 
most promising sources for a first direct detection of gravitational waves. Accordingly, a huge effort has been done 



00 . 

up to now in order to be ready in time for analysing the inspiraling binaries data delivered by the interferometric 
detectors : detection of the signal by correlation with suitable templates (matched filtering), see e.g [j^-fiol or by time- 

I < 

j^Tjl frequency analysis [ll|Jl^| , and estimation of astrophysical parameters (mainly the masses of the stars and their spins) 



[|1 3| — 11 5[| - This part of gravitational-wave data analysis, concerning inspiraling compact binaries, is now well advanced 
and rather well understood. Besides, periodic sources such as rotating neutron stars are maybe as interesting as the 



binary inspirals, since, despite the low expected gravitational waves amplitudes, these sources have the advantage of 
being permanent. A number of studies has been also dedicated to the analysis of periodic sources : see e.g fit] or |f6| ] 
and references therein. 

Historically, supernovae (SN) were the first envisaged gravitational-wave emitters and first resonant detectors have 
been designed to be sensitive around the typical frequencies expected in such bursts of gravitational radiation, around 
1 kHz. With the construction of intrinsically broadband inteferometric detectors, this kind of sources has not been 
studied as much as inspiraling binaries or pulsars. 
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Expected sources of burst gravitational waves are first collapses of massive stars to neutron stars (type II SN) or 
to black holes. Modern simulations of the former show a rather small efficiency of gravitational radiation emission 
]l8| : amplitudes typically less than 10 -22 are expected for SN at the distance of the Virgo cluster. Nevertheless, 
observation of very fast pulsars in the Galaxy (such as the one in the Guitar Nebula ^2|) may indicate that, at least 
in some cases, the collapse can be highly asymmetric and provides much higher gravitational wave strain amplitudes 
p3[ . Estimates of gravitational wave amplitudes from the collapse to a BH reach similar orders of magnitude as for 
previous type II SN @§. 

Another possible source of gravitational wave bursts occurs during the merging of two compact stars, at the very 
end of the binary inspiral. If the inspiral signal for binary neutrons stars is well understood up to the 2.5 Post- 
Newtonian order |p5f , we know only little about the signal waveform corresponding to the merging phase itself, since 
its computation requires in particular fully relativistic hydrodynamical codes, although some semi-classical attempts 



have already been performed, see e.g. p6| . Some recent estimates 27 give a maximum amplitude h ~ a few 10 -21 at 
10 Mpc within a frequency range of 1-2 kHz. This is just the order of magnitude of the noise level at these frequencies 
for interfcrometric detectors in their initial design; that leaves some hope for a future detection. Concerning BH 
binaries, an ambitious program called The Binary Black Hole Grand Challenge Alliance (28| is underway to handle 
the very difficult task of computing the waveform of the merging phase. 

Damped oscillations of excited BH's or NS's, like baby born NS's (just after the collapse), can also provide gravita- 
tional waves with detectable amplitudes ||. The corresponding waveforms are not really bursts like, they rather have 
some coherent structure (they look like typically a damped sine). However, their characteristic frequencies, above 
hundreds of Hz to tens of kHz, and their short damping times, make them belong to the category of signals of interest 



in this article. Note that the frequencies and damping times are exactly known for Kerr BH 29 and the detection 
of gravitational waves emitted by such a perturbed BH could provide a direct measurement of both its mass and its 



angular momentum |30|; of course, in this case, matched filtering, with damped sine templates, is the more suitable 
method. 

All these gravitational wave burst signals have the following features : short durations from milliseconds to seconds, 
frequencies from ~ 100 Hz to few kHz and a large range of waveforms. Filtering of such short signals in the ouput of 
interferometric detectors should therefore be as general or robust as possible, and designed with (almost) no a priori 
knowledge on the waveforms; this prescription of course forbids the optimal (Wiener) filtering as used for inspiraling 
binaries. Such general filtering methods are then necessarily "sub-optimal" , in the sense that they are less efficient 



2 



than the optimal filtering. In this article, we concentrate on the filtering of one detector's output, which is the first 
step in a detection process, the second one being the reconstruction of the gravitational wave signal from the filtered 
outputs of (at least) three detectors. The second step has been already studied in detail |H],[j2), while the first has 
attracted so far little attention in the litterature. Here, we study three filtering ideas for the detection of bursts 
in the data of interferometric gravitational-wave detectors, two of them being very general and the third one more 
specific. These methods are namely a "counting" method, where we count the number of bins which are larger than 
some threshold in a certain window, then a method based on the autocorrelation function of the detector data and 
finally a filtering based on the correlation of the data with a peak generic function. For each of them, we develop the 
statistical properties : link to gaussian statistics, number of false alarms, threshold definition ... In order to quantify 
the performances of such filters, we use as gravitational wave signals the SN catalogue from Zwerger and Miiller 
available on the web |33| , and compute, as a benchmark, for each of them, the maximal distance of detection obtained 
by the three filters; as a reference, we compare to the maximal detection distance calculated by optimal filtering. For 
this purpose, we will use, as a model for the detector noise, the minimum of the VIRGO sensitivity, which occurs 
precisely in the range of frequencies of interest. 

Of course, any burst filtering is unable to distinguish a non-stationary noise from a real gravitational wave event : 
such filterings will be sensitive to transient noise as well as to gravitational bursts. The goal of burst filtering in one 
detector is then mainly to act as a trigger and select interesting data streams in order to investigate coincidences with 
other detectors. It will also be useful to identify and study non-stationary noise in a single interferometric detector, 
ultimately providing vetos or cleaning procedures for "known" non-stationary noise sources. 

Finally, it is stressed that a general filtering approach, such as those proposed below, will be sensitive to unexpected 
sources and therefore may provide some insight into new physics. 

II. BURST FILTERING : SOME IDEAS 

Since we know little about the expected waveforms of burst gravitational-wave sources a robust filtering is required. 
Since such a filtering is wished to work as an on-line trigger, it should be fast. We study three of such filters in 
the following. Throughout the paper, we assume that the detector noise is white, stationary and gaussian with zero 
mean. For the numerical estimates, we chose the flat (amplitude) spectral density to be h n ~ 4 x 10 _23 /\/Hz and 
the sampling frequency f s ~ 20 kHz, so the standard deviation of the noise is a n = h n \fj~s ~6 X 10~ 21 ; we will note 
the sampling time t s = l// s . The value chosen for h n corresponds approximately to the minimum of the sensitivity 
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curve of the VIRGO detector [p4[ ; around this minimum, the sensitivity is rather flat, in the range ~ [200 Hz,lkHz], 
which is precisely the range of interest for the gravitational wave bursts we are interested in. This validates then our 
assumption of a white noise; otherwise, we can always assume that the detector output has been first whitened by a 
suitable filter p5j . 

A. Bin counting (BC) 

The principle of this first filtering method is quite simple. A data stream of length T being given (so containing 
N = T x f s data), we count the number of data (bins) whose value exceed a certain threshold, say s x <r n , in unit of 
the noise standard deviation. The method is illustrated in Fig.l. It follows the prescription about no preconceived 
idea about the waveforms to detect. In the absence of signal, the noise being gaussian, the probability that a data 
bin Xj is larger than s x a n is : 

P(\xi\ > sa n ) = 2 / e -^=dx. (2.1) 

J s V 2-7T 



It is then straightforward that N c , the number of bins above threshold, follows a binomial distribution and the 
probability that N c = n is 



P(N C = n) 



m r / . x m r , . N n N- 



ll 



erfc 



1 — erfc 



(2.2) 



where erfc is the complementary error function. Setting p = erfc(s/\/2), the mean of N c is then fi c — Np and its 



standard deviation is a c = \J Np (1 — p). It is well known that the normalised random variable N c — (N c — fi c )/^c 
behaves like a normalised gaussian variable, as soon as Np > 5 and N(l — p) > 5 [^6). These conditions will be 
fullfilled in every situation of interest, so we will consider now that the random variable N c is well approximated by 
a standard normal one. 

Two parameters are involved in this method : the window length T (or equivalently N) and the thresold s. We 
will discuss the window length later, with the other filtering methods. On the contrary, the choice of s is a specific 
issue for this approach. First, s should not be too large, as we expect low amplitude signals. Then, s should not be 
too small because then the filter would become very sensitive to the noise fluctuations with the drawback of a huge 
number of false alarms; as an example, if s = 2, -P(|a:i| > 2er„) ~ 4.6%, giving in average 46 'counts' in a window 
of N = 1000 sampled data. The optimal value for s is evaluated as follows : we compute the average distance of 
detection for the supernovae signals of the Zwerger and Muller catalogue as a function of s; the calculation of this 



distance is explained with full details in section [II. Of course, many realisations of the noise are generated and the 



results are averaged in order to reduce the influence of noise fluctutations. The result can be seen in Fig. 2. We choose 
accordingly s ~ 1.7 in all the following, but any value of s in the range [1.4,2.0] would also be reasonnable (giving a 
loss up to 10% in the average distance of detection with respect to the maximum). 

B. Norm filter (NF) 

This method has been initiated by the remark that white noise samplings are uncorrelated, while this is in general 
not true for a gravitational- wave signal. So the autocorrelation A x (t) = J x(t)x(t + r)dt of the detector output 
x(t) should reveal the presence of a correlated signal hidden in a uncorrelated noise. However, two problems arise : 
the computing time (depending on the window length) and the choice of the detection criterion. The norm of the 
autocorrelation function seems a good one, since it should probably distinguish between noise and noise plus signal. 
But the squared norm of the aucorrelation of a gaussian discrete variable is generally not a gaussian variable itself : 
it becomes gaussian only for very long data windows, as we have checked, and so it requires prohibitive calculation 
times. We could live with non-gaussian statistics however, but we first require simplicity for these preliminary studies. 
This is why we turn our attention to the maximum of the autocorrelation, which is nothing but the squared norm of 
the output x(t), since the autocorrelation of any function is maximal at zero. For the N sampled data in a window 
of size T = N/ f a , the output of this filter is then 

i=N 

A=Y / ^) 2 , (2-3) 
i=l 

where x(i) is the i th data in the window. When no signal is present, x(i) is pure noise and, under our assumption 
of gaussian noise, x(i) 2 , being the square of a gaussian random variable, is a Chi-square random variable with one 
degree of freedom. A is then the sum of N such Chi-square random variables with one degree of freedom, which is 
a Chi-square random random variable with N degrees of freedom. Its mean is [la = N and its standard deviation is 
a a — V%N- If N > 30, the random variable (related to the norm of the windowed detector output) 

A = V2A - V27V - 1 (2.4) 

is very well approximated by a standard normal variable |36[ . This fits the simplicity requirement and we will then 
use the output of A as a filter. The only parameter for this filter is the window size N; it will be discussed in the 
next section. 
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C. Correlation with single pulses (PC) 



Since many of the expected waveforms present one or several peaks, it seems judicious to use single pulses as filters. 
These pulses are modelised with truncated gaussian functions like : 



fr(t)=e X p(-^j, (2.5) 

with t lying in the range [— 3r, +3r], so that the function is truncated at about 1% of its maximum value. The only 
parameter for this set of pulse filters is the width r. The discrete correlation between the data x(i), i — 1,...,N and 
the pulse can be written as : 

N 

P(N, k) = J2 x(i + k)f T ([i - N/2]t s ) (2.6) 

i=l 

In the absence of a signal, the output of the filter P(N, k) is a gaussian random variable, as a sum of gaussian random 
variables weighted by the pulse fonction. The standard deviation of P is simply the square root of the sum of the 
squared weigths : 

o P = ^M[i-N/2]t s )*, (2.7) 

which can be recast as 

1 r m " t 
o* P *- Mtfdtc^-- (2-8) 

ts J — Nt s ''s 

We can then define a signal to noise ratio (SNR) as the filter output normalised by the standard deviation op : 
P = P/o-p, where P is the maximum of the function defined in Eq.(2.6). 

D. Practical implementation 

The two first methods BC and NF are very easy to implement in practice, as we can write simple recurrence 
relations for the calculation of the filter outputs in a given window, as function of the filter outputs in the previous 
window. For instance, for the BC, the output of the filter in a window of length N starting at the m th datum x(m) 
and ending at the datum x(m + N — 1) is N c (m). The next filter output N c (m+1) is obtained by moving the window 
by one bin, namely starting now at the datum x(m + 1) and ending at the datum x(m + N). The relation between 
N c (m) and N c (m + 1) can be cast as : 

N c (m + l)=N c (m) + e(m+l)-e(m), (2.9) 



where 0(i) = 1 if the datum x(i) is above threshold s x er„ and 0(i) = if is below. 

Similarity, for the NF, in the window of length N starting at the datum x(ra + 1), the norm of the data is 
A(m + 1) = X)£Lm+i x {i) 2 an d is simply related to A(m) by A(m + 1) = A(m) + x(m + N) 2 — x(m) 2 . This recurrence 
relation between A(m + 1) and A{m) allows a very fast calculation of the filter output A(m + 1) to be performed from 
the calculation of the previous filter output A(m). 

One advantage of this practical simplicity for both methods allows the computation of filter outputs with a window 
moving from bin to bin, which is not always possible (and anyway not necessary) in case of correlations with a 
predefined lattice of filters. 



Concerning the PC, we have first to built the lattice of filter, depending on only one parameter, the gaussian peak 
standard deviation r. The parameter space is the interval [r m ; n , r max ] . The distance between two successive filters 
of the lattice is noted At and the problem is to estimate Ar, which is a priori a function of r. The output of the 
correlation between a gaussian peak filter f T and a "signal" g is : 

</„,>=* Max,, Ulpnm, (2,10) 

where K is a constant. If g is itself a filter of the kind f T1 it is easy to show that the maximal correlation is obtained 
for t' — 0. Following y37| , we chose Ar such that 

< fr, fr > — < It+Ati It > 



< frjr > 



<e, (2.11) 



where e is the allowed loss in the signal to noise ratio. Expanding the ratio of Eq.2.11 to second order in Ar/r leads 
to the simple inequality 

<4e. (2.12) 



V r 

If the gaussian filter f T of width r belongs to the lattice, then the next one / t +At is then built from 

Ar = 2rVe. (2.13) 

Starting from the first filter of width r m i n , it is then easy to built the k th filter : its width is rfe = (1 + 2 v / e) fc_1 r m i n . 
The total number of templates in the lattice is finally the maximal integer n t such that (1 + 2y / e)™ t_1 r m i n < r max . 

For example, Fig. 3 shows the distribution of the 117 templates in the interval [1 ms,10ms] for a loss in SNR e = 10 -4 
(the one we will use in the next section); the choice of e — 1CP 2 reduces the number of filters to 13 for the same 



interval. We notice that we can allow for a very low loss in SNR and still obtain a reasonable number of templates; 
this is due both to the fact that we deal with an one-dimensional lattice space and to the 'smooth' dependence of n t 
on e. For instance, the same very low value of e = 1CP 4 and a parameter space extended to the (physically possible) 
interval [1 ms,l s] lead to a lattice of only 349 templates, which is easy to implement for an on-line processing. 

E. Threshold and false alarms 

Since the outputs of the three proposed filters behave like standard normal random variables (when no signal is 
present), it is convenient to define the same detection threshold for all, with, consequently, the same number of false 
alarms produced. As we expect weak signals, this threshold has to be low. On the other hand, we can deal with a 
large number of false alarms; these spurious events can be processed and discarded later when working in coincidence 
with other detectors. The relation between the detection threshold rj and the rate of false alarms rf a for a gaussian 
random variable : 

1 [°° f -x 2 \ 
2 x — = / exp — )dx = r fa . (2.14) 



Pto h V 2 

A false alarm rate rf a = 5 x 1CP 7 (~ 36 false alarms per hour for a sampling rate f s = 20 kHz) corresponds to a 
threshold r\ ~ 4.75, while a false alarm rate 10 times smaller corresponds to r\ ~ 5.20. For the results presented in the 
following), the chosen threshold is rj ~ 4.75. 

For the two first filtering methods, BC and NF, the situation is however not so simple because the outputs of the 
filters in two successive windows are in fact strongly correlated. For example, for the BC filtering, the filter outputs in 



two successive windows (starting respectively at the datum x(m) and x(m + 1) are related by Eq.2.9; it is clear that 
N c (m) and N c (m + 1) are the same or differ at most by ±1. So if the detection threshold is exceeded in a number 
of consecutive windows, there is in general only one "event" . This leads to redefine what is a detected event : an 
event is said to be detected in some time interval [mit s , m^ta] if the filter output 0(m) exceeds the threshold r\ for 
all values of m in the interval [mi, m,2\ and is less than r\ outside. This is equivalent to define a 'correlation length' 
(mi — mi)t s for the event to be detected. 

III. DETECTION OF SUPERNOVAE 

In this section, we use the catalogue of simulated gravitational-wave signals emitted during supernovae collapses 
and computed by Zwerger and Mullcr pof , to implement and test the three filterings described previously, in a realistic 
context. 



A. The catalogue 



The catalogue of Zwerger and Miiller 1 33 1 contains 78 gravitational- wave signals generated by axisymmetric core 
collapses. Note that, in particular due to axisymmetry, these are purely linearly polarized waves (/i x =0). Each of 
them corresponds to a particular set of parameters, essentially the initial distribution of angular momentum and the 
rotational energy of the star core, in the collapse models of Zwerger and Miiller. The signal total durations range 
from about 40 ms to a little more than 200 ms. The gravitational wave amplitudes of the stronger signals are of the 
order h(= h+) ~ a few 10~ 23 for a source located at lOMpc; that leaves little hope to detect such events with the 
first generation detectors. Concerning the shape of the waveforms (see Fig. 4), Zwerger and Miiller distinguish three 
different types of signals p0[ . Type I signals typically present a first peak (associated to the bounce)followed by a 
ringdown. Type II signals show a few (2-3) decreasing peaks, with a time lag between the first two of at least 10 ms. 
Type III signals exhibit no strong peak but fast (~ 1 kHz) oscillations after the bounce. The fact that type I and 
type II signals are characterised by strong peaks validates the choice of the filtering by correlation with generic peaks 
in order to detect such events. The 78 signal templates in the catalogue are not equally sampled, so we have first 
re-sampled them by interpolation at the desired sampling frequency (/ s = 20kHz in our examples). 

B. Optimal filtering and maximal distances of detection 

Since the 78 signal waveforms are known, we can explicitely derive the optimal SNR provided by the Wiener filter 
for each of them, and then compute the maximal distance of detection. We will then be able to build a benchmark 
for the different filters by comparing their results (detection distances) to the results of the Wiener filter. In all the 
following, we assume that the incoming waves have an optimal incidence with respect to the interferometric detector. 

Let's call h(t) one of the 78 signals and h(f) its Fourier transform. The optimal signal to noise ratio po is given by 

^ 2 i^ J (3 - 1} 

where is the one-sided noise power spectral density (hence the factor of 2). The noise is assumed to be gaussian 
and white with a standard deviation related to the constant spectral density Sh = by a n = h n x i/fs — V^h x fs- 
The SNR then becomes 

Pl= 2 ~k I \Hf)\ 2 df=^ I Ht)\ 2 dt. (3.2) 

As previously, a supernova signal is detected by the Wiener filter if p > t], where rj is the same detection threshold as 
defined above. Fig. 5 shows the maximal distance of detection for each of the 78 signals. The mean distance, averaged 



over all the signals, is about 26.1 kpc, which is of the order of the diameter of the Milky Way. A few signals can be 
detected at distances beyond 50 kpc, the distance of the Large Magellanic Cloud (LMC). The most energetic signals 
(number 77 and 78) can be detected at distances as high as 100-120 kpc. It is clear that this class of signals will be 
detected by the first generation inteferometric detectors only if the supernovae occur inside our Galaxy or in the very 
close neighbourhood. 

C. Detection by the three niters 

1. Window sizing 

The window size N for the first two filters is a compromise between the need for not too small number of bins, in 
order to garantee that the filter outputs are well approximated by gaussian random variables, and the rather short 
signals we are looking for. The first constraint is easily fulfilled for the NF method, where N must be larger than 
about 30. For the BC method, we must have A^erfc(s/v / 2) > 5 and N(l — erfc(s/v / 2)) > 5; with the optimal choice 
of s = 1.7, so that erfc(s/v / 2) — 0.91, the two inequalities give respectively N > 6 and N > 56. The constraints are 
not very severe and require that N must be between 80 and 150, as it has been checked by Monte-Carlo simulations 
(looking for the optimal window size giving the maximal mean distance of detection). 

For the pulse filter PC, the window size is automatically fixed by the definition interval of the gaussian peak 
[— 3t, 3t] for a filter of width t. The window size for the correlation with the filter f T is then TV = 6r/ s and depends 
then on r. 

2. Detection distance 

The efficiency of each filter is measured by the maximal distance of detection for each of the 78 gravitational 
wave signals of the catalogue. This distance is obtained by averaging over many noise realisations in a Monte- 
Carlo simulation. We present the results in two ways. Fig. 6 shows the number of detected signals as a function of 
the distance of detection for the three filters. Results of Fig. 6 combined with the results for the optimal filtering 
(Fig. 5) are reported in Fig. 7 in a normalised way. The histograms in Fig. 7 show the number of signals detected as 
a function of the reduced distance of detection for the three filters; the reduced distance of detection is simply the 
distance of detection divided by the maximal distance of detection computed with the optimal filter. The means 
of the distributions are respectively 0.22, 0.26 and 0.34 for the BC, NF and PC filters; these can be seen as rough 
estimators of the efficiency of the filters. The histograms in Fig. 7 give also an idea of the sensitivity of the filters 
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to the waveforms of the detected signals. We note that the histograms corresponding to the BC and NF filters are 
much more concentrated than the one corresponding to the PC filter; this is particularly impressive for the NF filter. 
This means that the two first filters are rather robust and their efficiency does not crucially depend on the details of 
the gravitational wave signals. At the contrary, the larger dispersion of the last histogram (PC) indicates that the 
response of the PC filter depends much more on the gravitational waveform. 

The global efficiency can be measured as the mean detection distance averaged over the 78 signals. The results 
are reported in Table 1. The efficiency of the BC and NF filters are about one third of the efficiency of the optimal 
filtering (max. efficiency), while the PC filter (for e = 10~ 4 ) has an efficiency about 58% of the maximal efficiency. 



Filter 


Optimal 


BC 


NF 


PC 


Average distance (kpc) 


26.1 


7.8 


9.3 


15.2 



Table 1 : Detection distance averaged over the 78 signals for the different filters. 

We note that none of the BC, NF and PC filters is efficient enough to cover all the Galaxy on average. Several 
signals however can be 'seen' anywhere from the Galaxy and even beyond; in particular the signals 77 and 78 can be 
detected up to the LMC by any of the three filters. 

In fact, concerning the PC filter, the mean distance of detection depends on e, the allowed loss in SNR, and 
consequently on the number of filters in the lattice. Fig. 8 shows the mean distance of detection as a function of the 
loss in SNR. As e decreases, the mean distance of detection becomes closer to the maximal value ( a little larger than 
15 kpc). We also notice that for high values of e (low number of templates), for instance above e = 1% (13 templates), 
the efficiency of the PC filter is still well larger than the NF/BC efficiency. 

3. Computation time considerations 

Apart form the criterion of simplicity of the filters, we require also that the data processing with these filters has 
to be fast enough to be implemented on-line as a trigger. 

An analysis of one day of data has been performed (in order to check the validity of the redefinition of an event 
-see section II. E- and the number of false alarms) on a DEC Alpha workstation, in about 14 mn with the BC filtering 
and in about 25 mn with the NF filtering. So there is about a factor 100 (resp. 57) between the data stream duration 
and the time needed to process them by the BC filter (resp. the NF filter). These two filters can then be used on-line 
without any problem. 

The PC filter is no more problematic due to implementation in the Fourier space and use of FFT's. Moreover 
the correlations don't need to be computed in successive windows. Indeed a template f s (gaussian of width s) has a 
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certain correlation length (v2s), so it is possible to compute the correlation with the template f s every T s , where T s 
is related to the allowed loss in filter efficiency. In practice, one order of magnitude on the calculation time can be 
gained with a loss of 1-2 % of efficiency. 

IV. CONCLUSION 

We have studied in this paper three filtering methods with the aim of building on-line triggers for the selection of 
burst-like events in the data flow of intcrfcrometric gravitational-wave detectors. Such a filtering needs first to be as 
general and as simple as possible; these two prescriptions are well fulfilled by the two first : BC and NF filterings. A 
third filter has then be studied in order to quantify what could be the gain with a more specific filtering, i.e. using 
some a priori information on the signals to detect (namely here the fact that the supernova signals contain at least 
one peak of short duration) . The catalogue of supernova gravitational- wave signals provided by Zwerger and Miiller 
has been used in order to built a benchmark for the different filters in a realistic way, and to compare each of them 
to the optimal (Wiener) filter. 

The results we find is first that a general filter, such as BC of NF, has an efficiency of about one third (0.30-0.36) 
of the efficiency of the optimal filter (maximal efficiency for a linear filter). The peak correlator PC (looking for the 
peaks in the Zwerger and Miiller signals) has a efficiency slightly larger (58% of the optimal filter efficiency). We note 
also that the general filters seem to be more robust (less sensitive to the waveform details) than the peak correlator, 
because of the smaller dispersion of their (normalized) responses to the 78 supernovae signals. 

Concerning a practical implementation, each of the filters can be implemented on line and can be used as a trigger 
in order to select events for off-line coincidences with other detectors. 

The results for the detection of the Zwerger and Miiller signals are not very optimistic. Moreover, as we have 
assumed optimal incidence for the gravitational waves, the detection distances should be divided by \/5 for sources 
randomly distributed over the sky S . So only a small fraction of the supernova events can be detected anywhere in 
the Galaxy or beyond. This is not a surprise, and, anyway, not the main point of this paper. 

We finally notice that the two general burst filters BC and NF we studied are in fact non-linear filters (they can 
not be reduced to a correlation with the detector output). This may encourage people to develop and study such a 
class of filters in the context of gravitational wave detection. In the near future, we plan to study other filters based 
on the autocorrelation function and on wavelets analysis. We plan also to study the effect of a more realistic noise on 
the response of this class of filters, in particular the effect of non-gaussianity and non-stationarity. We keep in mind 
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that non stationary noises can be treated as 'signals' and the burst filtering can help to identify them and finally 
understand the detector behavior. 
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FIG. 1. Principle of the Bin Counting filter. The filter select all bins that are above some level (3<r n on the example); here 
a signal, starting at time 0.5 s has been superposed to the noise. 

FIG. 2. Efficiency of the BC filter as a function of the threshold level s. The efficiency values have been normalised to the 
maximum value. 



FIG. 3. Distribution of the 117 templates for the peak correlator in the interval [1;10] ms with the assumption e = 10 



FIG. 4. Typical waveforms for the type I, II and III supernovae signals in the Zwerger and Miiller catalogue. 

FIG. 5. Detection distances calculated with the optimal filter for the 78 signals in the catalogue of Zwerger and Miiller. 
About 6 signals can be detected at distances as high as about 50 kpc (the LMC distance). 

FIG. 6. Histogram of the number of signals as a function of the maximal detection distance for each of the three filters (• : 
BC, o : NF and * : PC). 

FIG. 7. Histogram of the number of signals as a function of the normalised detection distance for each of the three filters. 
The normalised distance is the detection distance divided by the corresponding maximal detection distance computed by the 
optimal filtering. 

FIG. 8. Efficiency of the Peak Correlator PC versus the loss in SNR e. The efficiency is measured as the mean detection 
distance for the 78 signals in the Zwerger and Miiller catalogue. 
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